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Abstract 

We study spherically symmetric static empty space solutions in R + e/R model of f(R) gravity. 
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We show that the Schwarzschild metric is an exact solution of the resulted field equations and 



consequently there are general solutions which are perturbed Schwarzschild metric and viable for 



solar system. Our results for large scale contains a logarithmic term with a coefficient producing 

■ a repulsive gravity force which is in agreement with the positive acceleration of the universe. 
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I. INTRODUCTION 



The current accelerated expansion of the universe is one of the biggest surprises in cos- 
mology This acceleration may be due to unknown energy momentum components by 
the equation of the state p = — p or some other mechanism. On the other hand it may be 
justified by modification of general relativity. The idea to modification of general relativity 
is fruitful and economic with respect to several attempts made to give explanation for this 
accelerated expansion of the universe and shortcomings of standard general general relativ- 
ity 2-6]. Relaxing the hypothesis that gravitational lagrangian has to be a linear function 
of the Ricci curvature scalar, R, as in the Hilbert-Einstein formulation, one can take into 
account, an effective action where the gravitational Lagrangian is a generic f(R) function. 
One of the initiative f(R) models supposed to explain the positive acceleration of expanding 
universe has f(R) action as f(R) = R — p 4 /R [7|. After proposing the f(R) = R — p 4 /R 
model it appeared this model suffer several problems. In the metric formalism, initially 
Dolgov and Kawasaki discovered the violent instability in the matter sector lal. The anal- 
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ysis of this instability generalized to arbitrary f(R) models [9|, [10| and it was shown that 
an f(R) model is stable if d 2 f/dR 2 > and unstable if d 2 f/dR 2 < 0. Thus we can de- 
duce R — p 4 /R suffers the Dogalov-Kawasaki instability but this instability removes in the 
R + p 4 /R model, where /i 4 > 0. Furthermore one can see in R — p 4 /R model the cosmology 
is inconsistent with observation when non-relativistic matter is present. In fact there is no 



matter dominant era 
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121 ] . However a recent study shows the standard epoch of matter 



had claimed 



domination can be obtained in the R + /i 4 / R model 12| . The authors of 
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that the f(R) = R — ^ does not pass the solar system tests. Since the Schwarzschild metric 
is consistent with observations of the solar system, spherically symmetric solutions of a vi- 
able theory of f(R) gravity for the vacuum should have a satisfactory limit of Schwarzschild 
metric. As yet, the discussion on the weak field limit of this f(R) = R + P> 4 /R theory 
of gravity is done and there are a few papers which claim different results 13^4201] . In the 
recent works, authors have discussed f(R) — R ± p 4 /R theory of gravity for the case of 
weak field limit and obtained perturbation terms in the metric potential components which 



are independent 
value conditions 



rom the localized mass as source of inhomogeneous curvature and initial 
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221 ]. This matter independency appears because only a special solution 



is considered for the trace equation. In present work we show that the Ricci scalar depends 
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on the initial value conditions and the metric can be reduced to the space time of Mincowski 
as the Central body's mass tends to zero. By assuming small Ricci scalar, the present paper 
analysis includes both weak field limit and strong gravity regime of spherically symmetric 
solutions of f(R) = R + e/R gravity model at the metric approach for the vacuum. We 
have shown that the Schwarzschild metric is an exact solution of the resulted field equations 
and consequently there are general solutions which are perturbed Schwarzschild metric and 
viable for solar system. 



II. THEORY AND ARGUMENT 

In this section we investigate spherically symmetric static solutions of the vacuum field 
equations in R+ e/R model of f(R) gravity. The action for f(R) gravity is 

S = j d A x^(^ + ^, (1) 

where we have set K = 8ttG = 1 and C m is zero for the vacuum. Master Equations are 
obtained as 




whit 



g; = r; - s;r/2. (3) 

Contracting the field equation (J2]) we obtain 



R/3 + e/R= -D(e/R 2 ), (4) 

and it will be used to simplify equations ([2]). The analysis of spherically symmetric static 
solution can be carried out using Schwarzschild coordinate: 

dr 2 = A(r)dt 2 - B(r)- l dr 2 - rW. (5) 

Substituting the trace equation (jlj) into equations ([2]) give rise to the following equations: 
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R\ = B° (6b) 

K = K ( 6c ) 

where the last equation is trivial. 

At first, it is obvious that taking R = leads to R v = 0. Since, equations §5§ have exact 
solutions of the Schwarzschild for the metric components, i.e. 

A(r) = B(r) = 1 - — (7) 

where M is conventionally a constant positive real number. 

Now, we seek perturbed Schwarzschild solutions for equations (J6]). Hence we perturb R 
slightly from zero to the extent that 

R? , „ 

|i| <<C L (8) 

This can be done by adding perturbations to the initial value conditions of the Schwarzschild 
in some place, namely in r = tq. In fact, since the first derivative of B and the second 
derivative of A with respect to r appear in equations [6] we should have three arbitrary 
constant real numbers in a general solution. Consequently we can expect that equations 
([6]) have a solution which is a perturbation of the Schwarzschild metric in a spatial interval 
containing r . Approximating the right side of equation fl6a|) under condition (jSJ) we have: 



R r r = -R\ (9a) 
R\ = R 9 e (9b) 



Writing the Ricci tensor components in terms of A(r) and B(r), i.e. 



A'\ 2 AB + rB'A' BA" 



B' B'A' ( A'\ 2 BA" , . 

- BA' rB' + 2B-2 , , 

fl 2 = + — SS— ■ < 10C > 
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equations 9 can be separated as follows: 

B' = — - (11a) 
r 

A' 2 A' A 

A" = — (1 + B) + — (B - 1), (lib) 

2A 2rB y ' r 2 B K h y ! 

where ' denotes derivation with respect to r. The solution to equation (11 lap is 

B (r) = 1 - — (12) 

where M is a positive constant expected to be close to mass of the central body. As we 
mentioned above, answers for equations (jfJJ) will be perturbations to Schwarzschild metric 
by choosing initial value conditions slightly deviated from the Schwarzschilds'. If such a 
solution exists, we can regard A(r) as B(r) + a(r), where a(r) is very small in comparisons 
to B(r). By assuming 



\a(r)\ < \B(r)\ (13) 



equation lllbl can be linearized as 



.. 2M(M-r) (r 2 -5Mr + QM 2 ) , 

a" = — — J-a - : ~7j — -a' (14) 

r 2 (r-2M) 2 r(r - 2M) 2 ' V ; 

with the following solution 

a(r) = bB(r) +ja(r), (15) 
where b and 7 are constant real numbers and 



a(r)=B(r)Ln 2^ + 2y/r -2M - y/B(r). (16) 

Finally the metric components are presented as follows 

B(r) = 1 - — (17a) 
r 

A(r) = (1 + b)B(r) +7«(r) (17b) 

Since bB(r) is an exact solution of equation II lbl and leads to R=0, then only 7a (r) should 
be regarded as a perturbation to the Schwarzschild metric. Using the obtained answers R 
is written as 



As one can see when r — > 2M, B(r) tends to zero faster than a(r) and the absolute value 
of R increases rapidly. In this case as long as 

7 2 « B(r), (19) 

solution (117bl) is a good approximate answer to equation (lllbl) . But for the obtained answers 
to be valid, we should also check the assumed condition in equation [HJ By inserting the 
statement for the curvature scalar (i.e. Eq JT8|) into equation [8] it necessitate 

' < 4r 3 (r - 2M) . (20) 

If we set 7 ^ 0, equations [T9l and [201 predict that the metric field can not be a perturbation 
to the Schwarzschild on the horizon i.e. r = 2M. 

On the other hand, when r — > oo , a(r) tends to infinity logarithmically and it can't 
fulfill the assumption of equation ([13]) up to spatial infinity. In fact for r ^> 2M, 

A(r) «d + C 2 Ln(r), (21) 

is a pretty good approximate answer to equation lllbl with C\ and Ci being constant real 
numbers. For the latter case, curvature scalar R takes the form 

*M « (22) 

Since for r S> 2M the curvature scalar tends to zero rapidly, the obtained solutions do not 
violate the condition R 2 /\e\ <ti 1 and equations fTTTT) are valid. 

III. CONCLUSION 

We studied spherically symmetric static empty space solutions in R + e/R model of f{R) 
gravity. We have shown that: 

• The Ricci scalar depends on the initial value conditions and by setting them to zero the 
central body's mass vanishes and the metric reduces to the space time of Mincowski. 

• In the limit R = 0, solution of the model is Schwarzschild metric 

• In the ^ < 1 for the intermediate region, the solution of the model is perturbation 
to Schwarzschild metric and shows where r — > 2M, the solution of model can not be 
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a perturbation to the Schwarzschild metric. Because the requirement n <C 1 is not 
satisfied. 

• In the large scale limit we have B{r) — > 1 and A(r) take a logarithmic form with 
an arbitrary coefficient constant. Whereas the scale of the universe is not infinite, 
so the logarithm term is definite. As well as in the cosmological scale, gravity force 
depends on the CVs sign in which C2 < generates a repulsive gravity force which is 
in agreement with the positive acceleration of the universe. 
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